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ABSTRACT 


Since the collapse of the Soviet Union in 1990, the United States Navy finds 
itself searching for and defining new roles in the post-cold-war period. Design 
aids previously utilized to analyze the Soviet threat must be redesigned to 
assist in timely crisis response and the build-up of appropriate force levels for 
regional contingencies. Currently, no such model for crisis response exists for 
an operational commander and his staff that will help select the forces to 
respond. The Employment and Crisis Response Model, developed to assist in 
the contingency planning process, provides an indication of the total number 
of days required for the build-up of the desired force level on scene. The rapid 
calculation of this information allows decision makers to quickly analyze 


different response options when faced with a crisis. 
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DISCLAIMER 


The reader is cautioned that the computer program developed in this 


research may not have been exercised for all cases of interest. While every 


effort has been made, within the time available, to ensure that the program is 


free of computational and logical errors, it cannot be considered validated. 
Any application of this program without additional verification is at the risk 


of the user. 
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I. CHANGING TIMES AND NEW CHALLENGES 


As the world evolves, so too does the National Security Strategy of the 
United States. The rapidity of change and the multi-polar threat requires a 
decision maker to make critical and timely decisions during crises. A tool to 
assist him would be of great benefit to meet any contingency or operational 
requirement. Currently there is no such tool available, therefore the need for 
such an instrument is evident. The ideal tool is an easy to use computer 
model that can be quickly incorporated into existing computer architecture. 

The purpose of this thesis is to discuss the importance of a presence and 
crisis response model, and to provide the basic idea and programming logic to 
develop the model. A sample model, limited to aircraft carriers (CV) using 
hypothetical employment schedules is developed. Of particular interest is the 
determination of “optimal” employment strategies based on current Flexible 
Deterrent Options (FDOs) and available assets. Response options are based on 
calculation of response times, using a network shortest-path algorithm. 

This chapter discusses the need for a crisis response model and examines 
four different forward presence and crisis response models currently being 
utilized. Model limitations are also discussed. Chapter II introduces the 
Capabilities Assessment and Evaluation System (CASES) and includes the 
purpose behind the framework, current models being utilized in CASES and 
finally, planned additions to CASES. Chapter III discusses the basic idea 
behind the United States Navy Employment and Crisis Response Model 
(ECRM) including the complexities of presence, objectives of crisis response, 
the goals, measures of effectiveness, and finally, utility. Chapter IV is devoted 
to discussion of the program code, including specific information pertaining 


to the basic model. Finally, Chapter V provides an overview of the model and 


the model and concept, including limitations, extent of model testing, and 


areas of further study and future plans. 


A. THE WINDS OF CHANGE ' 

It is a time of great geo-political change. For the first time in the history 
of the United States, this nation stands alone as the only global superpower. 
The demise of the Soviet Union and Warsaw Pact in the early 1990s has, for 
the moment, decreased the likelihood of global warfare.! Although the 
nuclear threat continues to play a major role in political decision making, the 
dissolution of the Soviet Union has removed an enormous burden on the 
citizens and soldiers of the United States. Yet incidents in the Arabian Gulf in 
the 1980s and early 1990s are indicative of the relative instability of other 
regions of the world. Without the existence and direct influence of the Soviet 
Union, many nations or cultures may feel the need to address long 
suppressed external or internal objectives with military force, particularly 
during economically austere times. 

At home, the United States faces the daunting task of cutting the budget 
to reduce the deficit of the 1980s. The need to balance the budget is clear, and 
defense dollars and the “peace dividend” make convenient sources of 
additional funds for legislatures intent on showing progress in correcting 
budgetary woes. In this climate of decreasing budgets, the United States must 
develop new strategies to deal with diverse international crises such as the 
invasion of Kuwait, humanitarian missions in Somalia and Bangladesh, and 
the bloodshed in the Balkan peninsula. This is not a simple urdertaking, 
since the size of the Defense Department will be dramatically reduced. The 


future of naval strategy will hinge on flexibility and timely response to these 


1Sean O'Keefe, Frank B. Kelso, and Carl E. Mundy. “. . .From the Sea.” 
U.S. Naval Institute Proceedings, November, 1992, p 93. 


types of crises with limited assets. New computer models must be developed 
to assist planners in shaping strategy, preparing for unforeseen contingencies 


and maximizing force utility. 


B. COMPUTERS AND CONTINGENCIES 

Several computer models exist today to aid the planner in analyzing 
contingency operations, yet these models each address specific warfare 
mission areas (WMA). THOR provides data on strike warfare, Carrier Battle 
Force Logistics Expenditure and Resupply (CLEAR) models battle force 
logistical issues, and other models exist, or are being developed to address 
such WMAs as ASW (anti-submarine warfare), ASUW (anti-surface warfare) 
and amphibious warfare. Most models were originally conceived and written 
to deal with the Soviet threat. Clearly, such scenarios have lost their urgency, 
if not their importance. These models do have utility in the era of Major 
Regional Crises (MRCs), Lesser Regional Crises (LRCs) and Low Intensity 
Conflicts (LICs), however, one important aspect is ignored by these models. Of 
greater importance to the decision maker and analyst during a crisis are 
questions pertaining tc what forces are available to respond, and when those 
forces can arrive. A computer model which determines the location of 
available assets and calculates crisis response times allows alternative 
responses to be quickly generated and evaluated. Although computer models 
havo been developed for various studies, none are adequate for assisting the 
decision maker in this process. 

Several models that calculate forward presence and crisis response were 


examined for suitability. Two models were developed at the Center for Naval 


Analyses (CNA) to study forward presence and crisis response. Lieutenant 


2 Bruce Anderson and John Flynn. “CASES: A System for Assessing 
Naval Warfare Capability.” Joint Directors of Labs C3I Symposium, 1990. 


Commander Kevin Becker developed an APL (A Programming Language) 
model and John V. Hall developed a spreadsheet based model [Refs. 1 and 2]. 
Two other mode:is were examined; another spreadsheet model developed by 
Michael J. 1: .demann, written for a study conducted at the Naval Surface 
Warfare Center (NSWS), and one model developed for the Naval War 
College (NWC)/NSWC Naval Surface Warfare 2030 Symposium held at 
NWC in December 1991 [Refs. 3 and 4]. All four models are discussed briefly 


below. 


1. The Becker Model 

Becker developed a presence and crisis response model, written in 
APL, for a study conducted at CNA for the Commander-in-Chief, Pacific Fleet 
(CINCPACFLT). This model calculates the pres :nce and crisis response 
capabilities of Pacific-based aircraft carriers. Becker devised an algorithm, 
using the Pacific-based CVs, to examine their ten year deployment and 
maintenance cycles. It then calculates the expected presence of carriers in the 
Pacific for the next ten years. The resulting data is displayed as a histogram 
depicting the percentage of time various numbers of CV/CVBGs (aircraft 
carrier battle groups) may be present in a geographic region over time. 
Although Becker uses Pacific-based CVs for his analysis, any mix of surface 
platforms and submarines can be used in the model as long as tiie necessary 
historical employment data is available. Based on this presence data, he 
utilizes « transportation problem to determine estimated response times for 
various crises and CV peacetime deployment options. [Ref. 1] 

Since the model is primarily based on historical data, several 
weaknesses are inherent from the perspective of the theater level decision 
maker. Several assumptions, as stated by Becker, are necessary to determine 


the estimated presence of CVBGs. Two primary areas of difficulty are that CV 


schedules are assumed to be recurrent over the ten year period examined, and 
“make-ready-to-sail” times have been simplified.% 

Although CV schedules are based on some basic cycle plan, it is easy 
to understand that CV employment cycles are neither symmetric nor 
recurrent. Clearly, the last ten years of CV employments will not likely be an 
effective example for the next ten years.4 The CV has been the favorite tool for 
crisis response, and will continue to , lay a major role in diplomacy, however, 
anticipated reductions in the CV force will make it difficult to model future 
CV presence from historical CV data. In summary, Becker’s model cannot 
(nor was it intended to) be utilized as a tool for commanders in an 
operational environment. Another drawback for the decision maker is that 
the algorithm is written in APL, making modifications difficult for most 
users. The model’s strength lies in its ability to provide rapid “what-if” 


information for long-term general strategic planning and budget studies. 


2. The Hall Model 
Hall developed a desktop computer spreadsheet based presence and 
crisis response model at CNA. In it are even greater assumptions and 


simplifications.> The Hall model assumes fixed deployment schedules and 


3 Make-ready-to-sail time is the number of days required to make a CV 
that is in some type of maintenance availability deployable again. 


4 Several perturbations in the deployment patterns of CVs occurred in 
the last decade including responses to actions by Libya and Iraq. Although it 
can be argued that crises will arise in the future, it appears unlikely that CVs 
will be deployed again on the scale required to support Desert Shield and 
Desert Storm in the near future. 


5 While assigned to CINCPACFLT, CAPT William Schwabe, USNR, 
revised Hall’s model to reflect different assumptions and to model the 
problem more realistically. 


steady state conditions. Using average values, the model is able to calculate 
how much presence forces can generate in some geographic region, and how 
rapidly they could respond to crises. Again, this model provides data 
concerning average forward presence and average crisis response times, but is 
not intended to be utilized by the theater level commander in an operational 
environment in which position and readiness status are known at “go time.” 
It is intended to provide quick and easy answers to “what-if” questions 
concerning alternative deployment schemes such as the lengthening of 
deployments, alterations in maintenance schedules, forward basing of 
carriers, or changes in OPTEMPO ( operational tempo ) requirements, etc. 
Although useful to the Pentagon budget process and long-term strategic 
planners, neither serves the purpose of executing contingency operations. 
[Ref. 2] 


3. The Lindemann Model 

The Lindemann peacetime deployment model, similar in some 
respects to the Hall model, calculates an estimate of total surface force 
structure, by ship class, that is required to provide a specific worldwide 
presence. There are two versions of Lindemann’s model available, both of 
which were written for the Apple Macintosh™ computer using Microsoft 
Excel™. Version 1 calculates forward presence required based on geo-political 
issues. Since the primary purpose of this thesis involves crisis response, this 
version was not considered for evaluation. Version 2 allows the user to 
directly input presence requirements. This version also calculates various 
costs associated with the force level requirements necessary, and more 
importantly from the stand point of this thesis, also calculates expected CV 
surge capability. Unlike the previous two models, Lindemann does not limit 
his analysis to CVs, instead, distinguishing between CVBGs, ARGs 
(Amphibious Ready Groups), SAGs (Surface Action Groups) and CLGs 


(Combat Logistics Groups). He also analyzes ship types, such as CV, cruisers, 
destroyers, submarines and logistical support ships. Only his analysis of surge 
capability is limited to CVs. Nevertheless, this model suffers the same limits 
that the previous two models exhibit in that the results are expected values 


and averages. [Ref. 3] 


4. The Force Deployment and Estimation Model 
During the Naval Surface Warfare 2030 Symposium II, held at the 
Naval War College in Newport, Rhode Island, in December 1991, similar 
issues of presence and response were examined. Of particular interest is the 
Force Deployment Estimation Model (FDEM) developed for the symposium. 
FDEM allows users to: 
eExplore alternative deployment concepts and innovative force 
packages; 
eStudy force structure effectiveness; and 


¢Study trade-offs between alternate naval force deployment concepts.® 


Based on a specified force structure and world events, the model 
provides as output, a proposed allocation of forces. The inputs to FOEM can 
be adjusted to study hypothetical force structures, force packaging plans, 
basing configurations and demands for naval forces. As with the Hall and 
Lindemann models, however, this model represents a steady state response. 
FDEM addressed similar issues studied and analyzed by the Lindemann 
model, and may be the result of the efforts by NSWC to answer equivalent 
questions. [Ref. 4] 


6 Game Book. Naval Surface Warfare 2030 Symposium II, Newport, 3-5 
December 1991, p. V-1. 


C LIMITATIONS AND SOLUTIONS 

Each forward presence model discussed has strengths and applications, 
but none are designed for the theater level commander in analyzing options 
for specific contingency operations with specified requirements for force 
presence. The ECRM (Employment and Crisis Response Model) is developed 
in this thesis to assist the theater level commander. 

In ECRM, existing ship quarterly and cycle schedules should be used to 
provide real, timely and accurate data for the program. The drawback to this 
method is that conversion of employment schedules into a format suitable 
for computer calculation will be time consuming. Nevertheless, for 
contingency planning the best inputs must be used in providing options for 
planners and decision makers. 

In the current milieu of flexible crisis response, earlier stand-alone 
analytical warfare mission area (WMA) models, such as THOR, are too 
cumbersome to be of any use to the planner who must make critical and 
timely decisions. Other models such as the Becker and Hall models discussed 
earlier, are simply not applicable. The rapid development of computer 
technology, however, has enabled programmers to integrate individual 
WMaAs into one package. 

One such package under evaluation is the Capabilities Assessment, 
Simulation and Evaluation System (CASES) now in use at CINCPACFLT.’ 
There are, however, notable shortcomings; CASES models campaigns by 
assuming that assets are on station and available for operations. What if there 
are no assets available in theater? What are the available options to the 


decision maker? How soon can the planner expect to be able to conduct 


7Bruce Anderson and John Flynn. 


military operations in support of diplomatic and economic efforts? These 
crucial questions that are not answered by CASES. 

The importance of such questions are obvious in a multi-polar world. 
The draw down of military forces will lead to alternative deployment 
schemes by the Navy to meet global commitments. The planner must know 
the location and readiness of his forces , and when a crisis develops, he must 
be able to quickly determine the means of response with a preferred course of 
action. The need for a model to provide planners with just such options is 


evident in this remark by Dr. Ray Runyan, Senior Analyst at CINCPACFLT, 


A major initial requirement for timely and effective crisis response 
is the build-up of forces to appropriate levels, i.e., what forces are 
available, where are they, what is their status, when can they be on 
station, etc. CASES has much of this information now or has access, but 
a planning algorithm to bring it all together for the planner is needed.® 


To address the problem identified by Dr. Runyan, a model to assist in the 
planning phase of any contingency is to be developed. The United States 
Navy Employment and Crisis Response Model (ECRM) will show how to 
solve the problem posed and provide a valuable tool to assist the decision 
maker. Ultimately, the model can be integrated into the CASES framework to 


allow maximum flexibility and speed for the planner. 


8 Ray Runyan. “ FPC Models Board.” Models Review Board, San Diego, 
16 June 1992. 


II. CASES 


CASES originated as a battle management application in DARPA’s 
(Defense Advanced Research Projects Agency) Strategic Computing Program 
to demonstrate distributed heterogeneous processing, parallel processing, and 
knowledge-base support. It is a decision-support system used by CINCPACFLT 
to support the development analysis and execution of war plans and 


contingency operations.? 


A. THEBASICS 

In essence, CASES is a classic campaign analysis tool that is utilized in 
assessing warfighting capabilities on a theater level. In the past, campaign 
analysis was primarily conducted by operations analysts who modeled each 
separate warfare mission areas (WMA) as a batch computer process. Such an 
undertaking was formidable due to the tremendous amount of information 
that was required to be assembled, modeled in detail, and evaluated. Analysts 
were required to prepare model inputs, interpret model outputs, apply 
statistical principles and present their findings from multiple computer runs. 
Normally such a task required weeks or even months to accomplish. CASES 
automates the assessment process through an interactive workstation to 
perform analysis in hours or days.1° 

Standard Navy models are integrated within the CASES interactive 
framework which allows the user, typically military operators rather than 


analysts, to set up and visualize warfare plans and scenarios. This was 


9 Bruce Anderson and John Flynn. 


10 Thid. 


10 


accomplished by taking advantage of recent developments in GUI (graphical- 
user-interface) technology, as well as increased capabilities of workstations 
and local-area management systems and networks. Through a series of 
windows and selectable options, the user can access the huge database and 
create campaigns at the theater level, or simply analyze individual models 
within the CASES framework. With this capability, the operator can readily 
compare various contingency scenarios for the decision maker in minutes or 
hours rather than days or weeks. During the initial action and subsequent 
contingency planning phases, this capability is fairly clear. 
CASES release 6.2 is the latest model available and includes the 

following WMA models within its architecture: 

e ASW Environment (RAYMODE, PE, GEM-DB, etc.); 

e ASW Sensor Capabilities (Sweep/Search); 

¢Campaign-level ASW (NACM); 

eCVBF Defense (ASBAT); 

eAir Threat Assessment; 

eStrike (THOR); and 

¢CVBF Logistics Expenditures And Resupply (CLEAR). 


Additional models are planned for inclusion into the CASES 
environment. Those models are: 
¢War at sea (WASPS); 
¢Amphibious Warfare; 
Shallow Water ASW; and 


Land Campaign. 


B. CASES, FDOS AND CINCPACFLT 

CINCPACELT has begun to examine the concept of force packages as a 
means to simplify contingency planning and response. A classified message 
from CINCPACFLT contains detailed FDO (Flexible Deterrent Option) 
packages for various crises, devised from historical records involving force 
employment as well as diplomatic and economic actions taken by the United 
States in the past. The FDOs available to decision makers run the gamut from 
simple diplomatic actions, through various economic and humanitarian 
actions, to complete wartime force packages. Although general in nature, 
these FDOs are specific enough to allow the theater-level commander to 
deploy necessary forces as an initial response. 

Utilizing the FDO concept, it becomes apparent that combining FDOs 
with a crisis response mode] will allow the theater-level commander to 
examine available assets and select the best package to respond to the 
contingency. He will also have the ability to compare different FDO packages, 
if he chooses, by making simple user inputs and comparing the data. 

Enclosing this within the CASES framework increases the analytical 
capabilities further while maintaining simplicity. Not only will the decision 
maker be able to compare the FDO packages available and deploy the most 
suitable forces in the shortest amount of time, he will also be able to have 
CASES provide the necessary analytical data to examine, amongst other 
options, the strike, ASW or ASUW capabilities of the various FDOs 


suggested by the model, or input by the operator. 


III. US NAVY SHIP EMPLOYMENT AND CRISIS RESPONSE MODEL 
(ECRM) 


Four critical mission areas have been defined by ex-President Bush in 
shaping national strategy for the future [Ref. 5] and are echoed in the Secretary 
of the Navy’s White Paper, . . . From the Sea. Those four areas are strategic 
deterrence and defense, forward presence, crisis response and reconstitution.!! 
The requirement for a crisis response model is clear, and CASES provides an 
excellent framework to incorporate such a model. ECRM addresses the issues 
of forward presence and crisis response, making the model a valuable asset 


and tool for the planner and decision maker. 


A. THE MODEL 

The objective of ECRM is to provide the decision maker such as 
CINCPACFLT with a valuable tool in planning for contingency operations in 
his theater of operations. The model determines all assets available to the 
planner, and from within a group of eligible candidates, selects the best choice 
of ships to respond to the crisis. The response is dependent on the nature of 
the crisis, the level required, and the CINCs choice of FDOs. The model is 
intended to be implemented with all major U.S. Navy surface combatants in 
the Pacific, however, the scope of the thesis is limited to the employment and 
response of CVs so that the model can be easily followed. The concepts and 
ideas demonstrated here can be extended to include all major ship types such 
as amphibious warfare assets as well as cruiser/destroyer type ships. Figure 1 


provides a flow diagram of the basic operation of the model. 


11 Stan Weeks. “Crafting a New Martime Strategy.” U.S. Naval Institute 
Proceedings. Jan 1992. p. 32. 
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When a crisis develops, the model operates with the following data: 
eShip sch. tule database; 
*Region of interest; 
*Nature of the crisis; 
eSpeed of advance (SOA); 
*Possible port stops; and 
°FDOs. 


With this data, the model obtains a “snap-shot” of the current situation 
which serve as the initial conditions of the model.12 

Next, the model selects several possible deployment and response FDOs. 
As an example, suppose that a typhoon devastates the island of Maui. The 
CINCPACFLT user model determines that the crisis requires the disaster 
relief or humanitarian mission FDOs. From within this sub-group, the extent 
of the required assistance is determined, based on information provided by 
the user. An austere option could include sending a single destroyer to the 
island to provide medical and communications assistance as well as some 
limited fresh water and emergency power requirements. If the crisis is of 
greater severity, the next option would be sending a task unit to the area, 
centered on a major amphibious vessel with LCAC capability, as well as Navy 
construction battalion (CB) detachments. For this thesis, the required force 
will be the number of CVs required to be on station. Once the FDO is selected, 
the algorithm examines all available assets. Initially, the algorithm sorts the 


assets by status, i.e., whether the CV is conducting operations nearby, is 


12 For the purposes of this thesis, a hypothetical schedule based on ship 
homeport and specific deployment regions was utilized as the data for the 
model. It can be assumed, easily enough, that the data originated from 
schedules provided by the various fleet scheduling conferences. 
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currently in an extended maintenance availability or simply returning from a 
deployment. From within this sorted group, the algorithm selects the CV(s) 
with the shortest response time. The shortest response times include normal 
port stops scheduled for replenishment purposes and any necessary make- 
ready times. The initial output display is in tabular form so that the decision 
maker or aide can rapidly determine if the CVs selected are in fact, the best 
options. 

Once the desired FDO is selected, the model chooses the best assets to 
respond, calculates the required measures of effectiveness (MOE) and outputs 
the data to the user. The output data consists of the “optimal” FDO with the 
required MOEs. In addition, the model also provides data for FDOs that are 
one level above and below the recommended package to provide the decision 


maker with additional options. 


B. MEASURES OF EFFECTIVENESS AND MODEL OUTPUT 
There are two primary measures of effectiveness which are used to assess 
the crisis response requirements: 
*Time delay for the complete FDO to arrive at the scene. 
«Amount of combat potential or humanitarian assistance level 


available in the region each day during the buildup. 


The output of the program provides information to the user on the 
amount of combat potential/humanitarian assistance available as a ratio of 
numbers on station each day versus the user specified total combat potential 
or assistance level required to meet the crisis. With this information, the 
user can readily determine when all required units of the FDO will be on 


station, and observe the number of days required to reach this level. 
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C. AN EXAMPLE 

To provide an example on how the model runs, suppose a crisis 
develops in the Philippine Islands. A coup d’etat has toppled the 
democratically elected government. At the moment, CINCPACFLT has six 
CVs that are assigned to the Third and Seventh Fleets. CV1 is in an upkeep 
period in Alameda. CV2 is on station in the Indian Ocean, monitoring 
activities in the Arabian Gulf. CV3 is nearing the completion of an upkeep 
period in San Diego. CV4 is currently undergoing an overhaul in 
Washington and is not expected to be available for more than 30 days. CV5 is 
underway conducting training operations in Southern California. CV6 is 
underway conducting local operations near Japan. Based on the information 
obtained, the model outputs Table 1. 


TABLE 1: FDO ALLOCATION 


FDO Allocation 


Note: Estimated Days include standara !agistics stops and make-ready 


17 


From the table, CV6 appears to be the best option, since it is already at 
sea. Assuming an average transit speed of 18 knots, CV6 will require 3 transit 
days to arrive on station. Since the FDO only requires a single CV, CV6 would 
be the “best” choice. The next choice would be CV5, underway in Southern 
California. With a 14 day transit, plus one day for supply on-loads in Pear] 
Harbor, she could make it on station in 15 days. CV3 is nearing the end of an 
upkeep period with only one day remaining. With transit time and stops 
included, she would require 16 days. CV1 is also in an upkeep period, but 
would require more time to prepare for a deployment. CV2 is not permitted 
to leave the Indian Ocean, therefore, her estimated response is expected to be 
greater than 30 days. CV4 cannot even be considered for the response, since 
completion of her availability would require more than 30 days. 

At this point, the decision maker, CINCPACEFZT, has the option of 
changing the recommended order of response priority. Alihough the FDO 
package does nct require a second CV, she decides that another CV would be 
prudent. She decides that the model solution is not the best choice because of 
the excessive amount of time required for the second CV to respond to the 
crisis. Instead, she consults with her staff and decides that CV2 will be the 
second CV to respond. Orders are dispatched to CV2 ordering her to leave the 
Indian Ocean and proceed to the Philippines. 

The recommended priority is then adjusted to reflect the new FDO 
allocation solution, and the analysis begins. The modei calculates the 
percentage of the FDO to arrive on station from D-day and outputs the time 
delay required for the entire FDO to arrive on station. In addition, the model 
also outputs to an external text file, the cumulative percentage of the FDO 
arriving on station each day. This data can then be displayed graphically using 
a separate commercially available graphics generating computer program 


such as DeltaGraph Professional® or any good spreadsheet program. The 
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output of the example problem is provided as Figure 2, graphed with 
DeltaPoint’s software DeltaGraph Professional® [Ref. 6]. 

The model selects the next lower FDO, in this case, no CVs (ultimately 
this would imply a SAG, but one CV is the default) and repeats the 
calculations. Then, the model selects the next higher option, in this case, two 


CVs, and completes the analysis. 


| OU a 
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Number of Days 


-m— FDO(L) 


-e@— FDO 
—4— FDO(H) 


Figure 2: Graphical Display of Percentage of FDO Arriving on Station. 


IV. PROGRAM 


The goal was to write the computer code, included as Appendix A, ina 
manner that would be easily transportable to the UNIX environment, or 
even the PC desktop environment.!3 This necessitated the use of ANSI 
(American National Standards Institute) conforming C computer code. Since 
the C programming language is fairly prevalent and popular amongst current 
programmers, and the code conforms to ANSI C, the code can be easily 
modified for any future use for the model. This allows ECRM to be 
incorporated into existing theater level campaign analysis models used by the 


Navy or possibly into campaign analysis models used by sister services. 


A. DISCUSSION OF THE PROGRAM CODE 

Before discussion of the actual computer program can begin, a brief 
introduction to the data structures is provided. The main data structures 
utilized in the program are ship type, nodeSet type and edge type structures. 
The most important of these structures, the ship type, is used in the ship 
response priority and represents individual units, in this case, CVs. Later in 
this chapter, the description of the ship response priority, stored as a doubly- 
linked list, will explain the functions of the fields in the ship type structure. 


The ship type structure is made up of the following fields: 


eID unit identification; 


13 The reader is cautioned that the computer programs developed in this 
research may not have been exercised for all . a..-s of interest. While every 
effort has been made, within the time availatie ‘\ ensure that the programs 
are free of computational and logical errors, thy cannot be considered 
validated. Any application of these programs without additional verification 
is at the risk of the user. 
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*location location of unit; 

*status current unit status; 

ecP combat potential; 

ehv humanitarian value; 

edistance distance, in nautical miles, to the crisis; 
estops number of port stops normally required; 
edays number of days needed to respond; 
epNext pointer to next unit; and 

°pBack pointer to preceding unit. 


Although the combat potential and humanitarian value fields are 
available, those values were assumed to be uniform and equal (1.0 ) for each 
individual notional CV examined in this thesis. In reality, the values vary 
with ship class and type, and in some cases, individual ships within a class, 
depending on the capabilities inherent in specific units. In the future, 
quantification of both values will be required for calculation of the second 


MOE once all surface ships are implemented into the model. 


1. Entering Ship Locations and Status 
Following the initialization of variables and arrays, the main 
segment calls a sub-routine to read in the location data of assets in the theater. 
The external data file used is stored as an ASCII (American Standard Code for 
Information Exchange) text file in matrix form.!4 Each row of this matrix 
represents an individual unit and the columns represent time periods. Each 
time period is assumed to be two weeks in length, therefore, the entire matrix 


of unit locations represents approximately one quarter of the yearly schedule. 


14 Commercially available software such as spreadsheet programs can be 
used to enter or modify this data, and all other external file data, as necessary. 
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The program reads in each location, assigning a memory location for each 
entry. As each location is read in from the external file, the program builds a 
location array in memory similar to the external data matrix. 

Once the location information is entered, the computer then opens 
two files, one containing ship status files, and the other containing estimated 
make-ready times, in days. As the computer reads in the ship status code, 
memory large enough to hold the ship type data structure is allocated for each 
individual ship. When memory is allocated, the computer enters the 
individual unit status into the ship status field. The make-ready time is then 
determined for the individual unit and stored in the ship field representing 


days. A description of location and status codes is provided as Appendix B. 


2. Determining the Required FDO 
The program queries the user concerning the crisis code and source of 

the crisis. The crisis codes utilized in this model are provided as Appendix C. 
When the crisis code is entered by the user, the program calls the FDO 
function to determine the standard FDO required for the crisis. The integer 
value that is returned is used to determine the number of particular ship 
types required to respond to the crisis. For this thesis, the ten’s digit location 
indicates the number of CVs required, and the one’s digit indicates the 
number of LHXs ( LHA - amphibious assault ship, general-purpose; LHD - 
amphibious assault ship, multi-purpose ) required. Again, the program only 
examines the CVs required although the FDO code that is returned can be 
easily adjusted in the future to reflect other platform types, such as cruisers 
and destroyers. The program then queries the user as to the location of the 


crisis. 


3. Dial’s Algorithm 
Given the location and status of available assets and the crisis region 


of interest, rapid response times are a simple function of distance and speed, 
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as taught in basic algebra. If speed is fixed, then transit time is a linear 
function of distance, and the response time can be solved as a single-source, 
shortest-path network problem. It will be assumed that all assets will respond 
at a standard fixed SOA (speed of advance), assumed to be 18 knots. In the 
future, it may be necessary to incorporate variable speeds into the model if air 
assets are to be included. 

Based on the location of ships and the location of the crisis, the 
algorithm uses a network model to determine shortest-path distances 
between the regions of interest. The location of assets, the crisis region and 
way stops can be defined as nodes in a graph. The weights assigned to the 
edges of the graph are the geographic distances between the nodes, 
determined by DMA (Defense Mapping Agency) and provided in Table 2 
below. Using such distances, a graph can be created such as the one depicted in 
Figure 3 below. Although seemingly simple, the problem of shortest-paths 
rapidly increases in complexity with increasing numbers of nodes, which can 
be added in the future to reflect various geographic regions. It is anticipated 
that no more than 15 nodes would be required to provide sufficient 
resolution to encompass the Pacific theater of operations. 

Several algorithms such as Dijkstra, label-correcting, Bellman-Ford, 
and Ahuja-Mehlhorn-Orlin-Tarjan are available to solve the shortest-path 
algorithm. In selecting a method to solve the problem, it is clear that a 
method for solving a graph with nonnegative weights and no negative 
weight cycles would be a sufficient algorithm, since geographic distances are 
never nonnegative. Dial’s shortest-path algorithm, a modified version of 
Dijkstra‘s algorithm, was utilized to calculate the shortest-path distance with 
the crisis region serving as the source node. Dial’s algorithm, as with 
Dijkstra’s algorithm, requires that all edges of the graph have nonnegative 
weights, and the graph must not have negative-weight cycles. The main 


element of Dial’s algorithm utilizes a relaxation technique. In short, 
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relaxation is a method that repeatedly decreases an upper bound on the actual 


shortest-path weight of each vertex until the upper bound equals the shortest- 


path weight.15 


No comparison of the various methods was conducted. A run time 


analysis of the various algorithms was not investigated under the 


assumption that the number of nodes would not significantly impact the 


results. Though Dial’s algorithm was utilized here, any of the other methods 


available to solve shortest-path problems could be implemented and used. 


FIGURE 3: A Simple Transportation Network 


Based on the distance solution obtained from Dial, the program must 
update the distance field for each unit. For each ship, ECRM compares the 
ship’s present location to the location of the first node in the distance 
sequence. If the locations match, the ship data structure is updated by entering 


the distance in the ship distance field, and calculating the number of days 


15 Thomas H. Corman, Charles E. Leiserson, and Ronald L. Rivest, 
Introduction to Algorithms. (Cambridge: The MIT Press, 1990), p. 520. 


required to respond to the crisis. Response days are dependent on the ship’s 
status. If the status code indicates that the ship is unavailable to respond, the 
number of days is set to thirty, the default for units not able to respond. If the 
ship is able to respond, the number of transit days is calculated by dividing the 
total distance by 24 (hours per day) and again by the crisis transit speed of 18 
knots. The number of days required for logistical stops while transiting to the 
crisis is added to produce the total steaming days required to respond. This 
value is then added to make-ready time to arrive at the estimated total 


number of days required to respond. 


TABLE 2: DISTANCES BETWEEN NODES 
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Defense Mapping Agency chart, The World , edition 1-DMATC, stock no. 


SF - San Francisco/West Coast of United States 
HI - Hawaii/Mid-Pacific 


YK - Yokosuka, Japan/Northwest Pacific 
PI - Philippine Islands/Southwest Pacific 
SI - Singapore 

IO - Indian Ocean 
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4. Sorting the Ship Priority List 

Following the calculation of total response times for each ship, the 
ships are then organized into a ship priority linked-list. Before the ship 
priority list can be displayed for the user, it must be sorted by status and 
distance. The ship list is reorganized as a binary tree array where the first node 
in the array represents the ship with the highest status value. This value is 
then exchanged with the last ship in the list, and the heap sort algorithm is 
called, recursively if necessary, to sort the new list. Again the first node is 
exchanged with the last node, and the process is repeated until all units have 
been sorted. The program then compares the status values in the ship list. 
Ships with equal status values are then sorted linearly to prioritize the listing 


by status and distance. 


5. Generating Response Table and MOEs 

The ship priority list is then output in tabular form, as in Table 1, so 
that the user can examine the response recommendations and make changes 
as necessary. The user has the ability to change the order of the response, the 
status of any ship, the number of days in port for logistical purposes, the 
make-ready time and the default transit speed. When changes have been 
entered, the ship priority list and table is updated, reflecting the new 
information. When satisfied with the resulting list, the user can resume the 


model, allowing for calculation of the measures of effectiveness. 
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V. CONCLUSIONS 


Recent upheaval of the global geo-political climate, combined with the 
capricious nature of the Earth, regularly demonstrate that crises around the 
world can occur suddenly and unexpectedly. In an era of diminishing budgets 
and downsizing, the need to respond to crises has not abated. In fact, it can be 
argued that with the United States as the sole superpower, our global 
responsibilities have increased. The United States Navy must be ready to 
meet these demands quickly and efficiently to ensure the restoration of 
security and peace. 

To assist the decision maker in determining an adequate response to 
crises, a decision aid that can analyze various FDOs and deployments with a 
campaign analysis model such as CASES is invaluable. Currently, no such 
deployment model exists to assist the planner and decision maker. 

Four presence and crisis response models were evaluated and 
determined to be unsuitable for analyzing actual contingencies. Therefore, a 
need exists for a contingency planning model that can be incorporated into 
the CASES framework. The developed Employment and Crisis Response 
Model is the solution. ECRM allows computer operators and watchstanders, 
rather than the analysts, to quickly determine courses of action in crisis 
response. 

One of the model’s strengths is its flexibility, allowing operators to 
examine the best available forces to meet a required FDO recommended by 
the model , based on Operational Support System (OSS) data for current 
locations of CVs and normal transit routes. If the model solution is not 
satisfactory, the user can alter the model solution prior to calculation of MOEs 
to meet the needs and desires of the decision maker in choosing the “best” 


force. This flexibility allows users to compare different deployment options 
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quickly and provide prompt answers to questions concerning alternative 


deployment options. 

ECRM is an adaptable model, since the program was written in ANSI C. 
This language also insures compatibility with CASES in the CINCPACFLT 
OSS. In addition, the computer code follows a modular design, allowing for 
easy enhancement of the program. Other service assets, such as Army 
divisions or Air Force squadrons can be easily programmed and added as a 
module by following the design of this thesis prototype. Although intended 
to serve as an input into CASES scenario analysis, ECRM also has the ability 
to stand alone on any decktop « ¢ notebook computer. 

With appropriate a. . from OSS and some modifications to the data 
entry interface, the CV model can be implemented at CINCPACFLT FPC as is. 


A. LIMITATIONS OF ECRM 

Although ECRM is definitely a push in the right direction, it is still a 
long way from providing the necessary capabilities essential in analyzing 
scenarios. Currently, ECRM is limited to calculating the response times of 
aircraft carriers. Useful as that may be, clearly it does not provide the 
necessary flexibility in describing full response capabilities, including FDOs 
made up of amphibious platforms, or surface action groups. 

ECRM does not determine the preferred FDO in a crisis. Simply put, the 
model does not examine and compare the various combat and humanitarian 
potentials of different ships and aircraft combinations to respond to crises. 
Doing so would probably necessitate the use of an integer programming 
model to obtain a solution. This does not, however, limit the usefulness of 
ECRM. The FDO concept was developed to allow decision makers to “plug- 
in” deployment modules that meet specific criuis criteria. These deployment 


packages greatly simplify the task of the decision maker, providing an 
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excellent basis for more detailed analysis of additional options. Once the FDO 
is selected, the model recommends the resources that will meet the 
requirements of the FDO in the shortest amount of time, based on the data 
and actual states of CINCPACFLT units. 

ECRM does not have the ability to process more than one crisis at a time. 
Nevertheless, the model has the flexibility to compensate. In most cases, the 
discrete nature of crises will allow allocation of forces for one crisis, followed 
by allocation to another. Analysis of one crisis will update the status markers 
and database, which will be reflected in the analyses of subsequent crises. 
Prioritizing crises prior to model input can greatly enhance the quality of the 
model output, as is expected. 

ECRM is not a substitute for thorough analysis conducted by expert 
analysts—their value in providing studies and modelling cannot be replaced 
by ECRM. ECRM only serves as a catalyst for rapid decision making when 
faced with a crisis requiring instantaneous action. 

In summary, the two main limitations to ECRM include the inability to 
determine “optimal” FDOs in a crisis, and the inability to handle 
simultaneous crises, however, the model does provide the necessary 


flexibility to allow users to minimize the effects of these limitations. 


B. FUTURE PLANS AND AREAS OF FURTHER STUDY 

This thesis does not develop an all encompassing CASES module that 
can be easily integrated into the framework immediately. The purpose of this 
thesis was to develop the basic model prototype required to program such an 
ultimate model, employing the ideas presented here. The goal of the model, 
as the title suggests, is to incorporate all major U.S. Navy surface combatants 
available to the decision maker and planner. By contrast, the prolific use of 


mines as a cheap substitute for capital ships has been demonstrated on 


29 


numerous occasions. It may be wise to include mine assets in a crisis response 
model. 

The eventual inclusion of combat logistics assets is recommended. 
Timely response of combat potential or humanitarian assistance is vital in a 
crisis, but immediate response is only a portion of the whole effort. The 
sustainability of FDOs in a crisis region with fuel, ammunition and supplies 
is of extreme importance both to the planner and decision maker. An integer 
programming algorithm or other method could be used to calculate the 
logistical requirements of each FDO, and then select the logistical assets 
required to sustain the FDO for some period of time is an important next step. 
CASES will have the ability to evaluate logistical concerns once the CLEAR 


(CVBF Logistics Expenditures and Resupply) model is fully integrated into the 
CASES environment. A means to tap into this model should be investigated. 

Modules for other Navy assets such as submarines, aircraft squadrons 
and construction battalions should be developed to enhance the capability of 
the model. To meet joint operational requirements, ECRM should include 
major Air Force and Army assets, such as light infantry or airborne divisions 
and rapid response tactical air lift. The drawdown in the Department of 
Defense will only increase the need for joint operations in crisis response. 
Modules for the sister services will be beneficial to the joint task force 
commander as well as the theater level decision maker. 

Finally, ECRM should incorporate a forward presence algorithm. Such 
an algorithm will allow planners and decision makers to assess the impact 
that certain deployment packages will have on future forward presence and 


short-term deployments. This “what-if” forecasting capability may provide 


additional constraints on selection of assets that make up the recommended 
FDO. 


APPENDIX A 


[PRONERER SERENE AH ERED EE TEER EEE EERE ESTE EREEREEEEEER NEES RENE OTE S CEES 


Structs Header Module 
Written by LT Takashi R. Yamamoto 
Operations Research Department, Naval Postgraduate School 


*- * * 8 8 


assignments 
ED EH HH Ob OF HE HE HE OF % HH BE HE OH 9 DE ED Ib Ob Ob Ob ob ED 5b DD HH HH DD 


» 
* 
» 
* 
Contains all struct types defined, as well as #define * 
* 
/ 


#define maxNodes 6 

#define numCV 6 

#define speed 18 

#define infinity INT_MAX 
#define maxEdge 4268 
struct Posit 


{ 
short period; 
short location; 
struct Posit *next; 


struct ship 


short ID; 

short location; 

short | status; 

float cP; /* measure of combat potential */ 
float hV; /* measure of humanitarian value */ 
short distance; 

short stops; 

short days; 

struct ship *pNext; 

struct ship “*pBack; 


struct nodeSet 


( 
short node; 
short distance; 
struct nodeSet *pNext; 
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struct 
short 


pNode 


struct 


edge 
short 
short 


long 
struct 


node 


struct 


nodeSet *pBack; 
bucket; 


nodeSet *pTr; 
nodeOut; 

nodeln; 

distance; 

edge “*pNext; 


edge “*pAdjNodes; 
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LEPESANT E SER TORS EEE E EERO E SA TRS SR OE AON EO NESS I EN SETAE SOT OS eee eS 


7 Global Variable File * 
* Written by LT Takashi R. Yamamoto “i 
* Operations Research Department, Naval Postgraduate School * 
* 


* 


Global declarations for ECRM. * 


SRETRERER SEEDER EEE EEE EERE EERE SERE EEL EREES OL EON ETNA RE HEN ERE DAS O / 


#include "structs.h" 


extern FILE *data; 

extern short source, D, i, crisis, region, xsitSpeed; 
extern short fDO; 

extern Boolean good; 

extern struct node  adjList[ maxNodes ]; 

extern struct pNode B[ maxEdge + 1]; 

extern struct nodeSet *P[ maxNodes + 1]; 
extern struct nodeSet *pFirst, *pLast, *setS, *sLast, *sMarker; 
extern struct nodeSet *pVertex; 

extern struct Posit CVLocation{ numCV J; 
externstruct ship *CV[numCV]; 


JPTHREhs SERENE EERERERERET EEE UEE EEE SESE NASH SER SEON SEAS E ESSERE SOE eS 


: Crisis Code File 
Written by LT Takashi R. Yamamoto 
Operations Research Department, Naval Postgraduate School 
Procedures called from Procs and Funcs. 


be modified to reflect user applicatoins 


POSSE OEOEHERADHAAFEEESEDEEEDEAESESPEESEHEES SORES OREREETEDEEE EEE DDH 


* 
» 
* 
Prints hypothetical crises onto screen. This module should 2 
/ 


void PrintCodes() 
/* 
Called from procedure DetermineCrisis. Prints various 
crisis options onto console window. 


BY A 


printf( "[ 1] \t Coup d'Etat \n" ); 

printf( "{ 2] \t Bombing of US/ Allied Embassy or Base \n” ); 
printf( “{ 3 ] \t Revolution in country of interest \n" ); 

printf( "[ 4 ] \t Civil War in country of interest \n" ); 

printf( "[ 5 ] \t Territorial Dispute of interest \n" ); 

printf( "{ 6 ] \t Hostage Crisis \n" ); 

printf( "[ 7] \t Attack on US forces \n" ); 

printf( “[ 8 ] \t Humanitarian Aid \n" ); 

printf("[ 9 ] \t Non-combatant Evacuation \n" ); 

printf( "{ 10 ] \t UN Maritime Interdiction \n" ); 


[SEER EEEEEER REESE EE REE EEE EERE ER EEE EER ELEEL EDEL SEDELEEE SEEDS ODS 
Presence Calculation Module 


Written by LT Takashi R. Yamamoto 
Operations Research Department, Naval Postgraduate School 


x #* + § B&B * 


* 
* 
Called from main to calculate current ship presence . 
* 
'é 


‘tH 3b 1 Ob Ob Ob HD HD OOF 1 Ob EO EE OE 


#include <stdio.h> 
#include <stdlib.h> 
#include "structs.h" 


extern struct Posit CVLocation{]; 
extern struct ship *CVII; 
extern short crisis, region; 


struct Posit “*LastPtr, *currentPtr; 


void BuildArray( struct Posit *currentPtr, short cv, short p, short k ) 
/* 


Called from procedure ReadLocation. Creates a matrix 
of surface ship location, based on the data in the 
external text file “location.” 


*/ 
(*currentPtr).period = p; 
(*currentPtr).location = k; 
(*currentPtr).next = NULL; 
if (CVLocation[ cv ].period == 0 ) 


CVLocation[ cv ] = *currentPtr; 
LastPtr = currentPt; 


( *LastPtr ).next = currentPtr; 
LastPtr = currentPtr; 
1 /* end ‘if loop */ 


} /* end , cocedure BuildArray */ 


void FlushFile( FILE *fp ) 
/* 


Called by various procedures. Flushes data stream. 
“f 
short ¢c; 


while( ( (c = fgetc( fp ) ) != '\n') && (c != EOF ) ) 


‘, 


} /* end procedure FlushFile */ 


void UpdateStatus() 
/* 
Called from procedure ReadLocation. Opens the external 
files containing status and make-ready data for each 
ship. Reads the appropriate data into the corresponding 
fields of the ship type data structure. 


*/ 


FILE “fp, *fp2; 
short _ state,cv, mrt; 


fp = fopen( "Status", "r" ); 
fp2 = fopen( "Make-ready", "r" ); 


if ( fp == NULL) 
{ 


printf( "The ship status file cannot be found. \n" ); 
printf( "Please ensure that the file ‘Status’ is in the current directory.” ); 
exit(0); 

); 


if (fp2 == NULL ) 

{ 
printf( "The ship status or Make-ready file cannot be found. \n" ); 
printf( “Please ensure that the file 'Make-ready' is in the curreni directory.” ); 
exit(0); 

); 


for( cv = 0; cv <= numCV - 1; cv++ ) 
{ 
fscanf( fp, "%d", &state ); 
FlushFile( fp ); 
fscanf( fp2, "%d", &mrt ); 
FlushFile( fp2 ); 
CVI cv | = malloc( sizeof( struct ship ) ); 
(*CV{[ cv] ).1ID=cv +1; 
(*CV[ cv ] ).status = state; 
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(*CVI cv } ).days = mrt; 
); 


fclose( fp ); 
fclose( fp2 ); 


/*“** read in Location array ***** / 


void ReadLocation() 


/* 
Called from the main program. This procedure opens a text 
file containing the locations of ships within the area of 
interest. The completed matrix will have the presence of 
ships the current period, and for the next eleven periods 
( approximately one quarter ) where each period is two 
weeks in length. 
*/ 
{ 

FILE “fp; 

short k,cv,p; 

char ¢; 


fp = fopen( "Location", as on ); 
if (fp == NULL ) 
{ 
printf( "The ship location file cannot be found. Please ensure \n" ); 


printf( “that the file ‘Location’ is in the current directory.” ); 
exit(0); 


_— 
= 


p=0; 
cv =0; 
while ( ( c = fgetc(fp)) != EOF ) 
{ 
while ((c != '\n' ) && ( ¢ != EOF )) 
{ 
if((c!='Nt') && (ct="')) 
{ 
++; 
k = ChrToDigit( c ); 
currentPtr = malloc( sizeof( struct Posit ) ); 


BuildArray( currentPtr, cv, p, k ); 


); 
c = fgetc(fp); 


CV++; 
p=0; 


fclose( fp ); 


UpdateStatus(); 
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f TEA TESTER SESS RARE AES E ESAT TESTE RES REST SER EERE SEARS SSR ESE SS St 
* FDO Determination Module : 
" Written by LT Takashi R. Yamamoto * 
* Operations Research Department, Naval Postgraduate School * 
* » 
* Based on user entered crisis code, determines the i 
5g generic FDO package required to meet contingency. > 


FEROAEEREAREES ERE SEEELELEEEEYESEEEEREEE ERS ESEEREEAEHEEEEEREERERR / 


#include <stdio.h> 
#include <stdlib.h> 
#include "Globals.c" 


short FDORequired( short criCode ) 
/* 

This function will return the FDO code, based on the crisis 
code that is passed to the function. The FDOs used in this 
program are HYPOTHETICAL ONLY and do not reflect the actual 
codes utilized by CINCPACFLT. This is insure that the program 
remains unclassified. For realism, programmers should alter 
the function as required. of f 


short level, code; 


switch( criCode ) 
{ 

case 5: 

- Territorial Dispute with one friendly involved. 
Obviously, the response to this situation would 
vary with the level of force deemed necessary and 
the nations involved. For the purposes of this 
project, it will require 2 CVBGs. */ 


code = 21; /*2CVBGs, 1 ARGs */ 
break; 


case 6: 

/* US citizens in foreign nation taken as hostages by host nation 
Initially, this situation would require a 2 CVBG and possibly 
and ARG, however, for the purposes of this project, only the 
CVBG force will be deployed. */ 


code = 21; /* 2 CVBGs, 1 ARG */ 
break; 


case 7: 
/* Military attacks on US forces overseas */ 


code = 32; /*3CVBGs,2 ARGs */ 
break; 
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case 8: 
/* Natural Disaster Relief/Humanitarian Aide */ 
printf( "What level of relief is required?: " ); 
scanf( "%d", &level ); 
switch( level ) 
{ 
case 1: 
/* Level One - the lowest */ 
code = 11; 
break; 
case 2: 
/* Level Two - mid-level */ 
code = 12; 
break; 
case 3: 
/* Level Three - catastrophic */ 
code = 13; 
break; 
default: 
code = 1; 
}; 
break; 
case 9: 
/* Non-combatant Evacuation */ 
code = 11; 
break; 
case 10: 
/* Enforcing UN trade sanctions on SLOCs */ 
code = 30; 
break; 
default: 
ial The default value represents a crisis that requires only a 


a single CVBG. */ 


code = 10; 
); 


return( code ); 


[RHSSOAERRECEREENE SER OREE ER EES ES SS SHON ES TESOL NN ESEN OSE R RENTS ESS 


» Procedures and Functions Module * 
* Written by LT Takashi R. Yamamoto . 
* Operations Research Department, Naval Postgraduate School * 
* » 
i General procedures and functions used by ECRM. : 
* Except some functions for Dial's Algorithm, which remain * 
7 in the Dial module * 


HESEPREERARENER ESSE RERRASRLAORE NET EN RSE EERE EERE LASER ESORNS OER RO: /’ 


#include <stdio.h> 
#include <stdlib.h> 
#include "Globals.c" 


short portStops{ maxNodes ][ maxNodes J; 
struct ship “shipQ, *pQIndex; 


[RERPOREREE READ NEARER ERED / 


void Initialize() 
{ 


short i; 


for( i = 0; i <= maxNodes; i++ ) 
adjList[ i ].pAdjNodes = NULL; 
for( i = 0; i <= maxEdge; i++ ) 
B( i J.pTr = NULL; 
setS = NULL; 
good = false; 
xsitSpeed = speed; 


} /*** end of Initialize ***/ 


short ChrToDigit( short khr ) 
/* 
Converts character types to integer types. 


return( khr - ‘0' ); 


void DetermineCrisis() 
/* Although this program requires the user to manually enter 
the crisis code and region code, the final program should 
also be able to obtain this data from CASES directly. 
af 


printf( "US Navy Ship Employment and Crisis Response Model \n\n\n\n" ); 
PrintCodes(); 
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void 


while( !good ) 
{ 


printf( "\nEnter crisis category. " ); 

scanf( "%d", &crisis ); 

good = ( crisis >= 0 ) && ( crisis < 11 ); 

if( !good ) 

printf( “Invalid entry. Try again. \n" ); 

); 
printf( "\n" ); 
good = false; 
fDO = FDORequired( crisis ); 


while( {good ) 
{ 
printf( "Enter crisis region: " ); 
scanf( "%d", &source ); 
good = (source > 0) && (source <= maxNodes ); 
if( 'good ) 
printf( "Invalid entry. Try again. \n" ); 
); 
printf( "\n" ); 


DetermineStops() 
* 


Called from procedure DetermineDistance. Inputs the 
port stops between different nodes and stores the 
data in a two-dimensional array. 


*/ 


FILE “fp; 
char ¢ 
short k,x=0,y=0; 


fp = fopen( "Stops", "r" ); 

if (fp == NULL ) 

{ 
printf( "The port stop file cannot be found. Please ensure \n" ); 
printf( "that the file ‘Stops' is in the current directory.” ); 
exit(0); 

}; 


while ( (c = fgetc(fp)) != EOF ) 

( 
while ((c != '\n' ) && (¢ != EOF )) 
{ 


if (cc !='\t) && (c!='')) 


{ 
k = ChrToDigit( c ); 


41 


portStops| x I y ] = k; 
y++; 
}; 
c = fgetc(fp); 
} 
X++; 
y=0; 


); 
fclose( fp ); 


void DetermineDistance() 
/* 
Called from main. Updates ship fields, based on the 
distances calculated by Dial. Updats the location, 
status, days and stops fields for each ship. 
a 


struct nodeSet *pDistanceQ; 
short c, posit, steamingDays; 
Boolean found; 


DetermineStops(); 
for(c = 0;c <= numCV - 1; c++) 
{ 
pDistanceQ = setS; 
found = false; 
posit = CVLocation{ c ].location; 
while( !found ) /* no cv in region */ 
{ 
if( pDistanceQ == NULL ) 
{ 


printf( "CV not in this theater. \n" ); 
(*CV[ c] ).status = 4; 
(*CV[ c] ).days = 30; 
found = true; 
) /* If unable to locate CV in theater */ 
else 


found = ( ( *pDistanceQ ).node == posit ); 
if( found ) 
( 
(*CVI ¢ } ).location = posit; 
(*CVI ¢ ] ).distance = ( *pDistanceQ ).distance; 
steamingDays = ( ( ( *CV[ c ] ).distance )/xsitSpeed ) / 
24; 


if(((*CV[ cj ).status == 4 ) && ( posit != source ) ) 
(*CV[ c] ).days = 30; 


else 


(*CV[ c] ).days = (*CVI c ] ).days + 


steamingDays; 
(*CV[c ] ).stops = portStops| source - 1 }[ posit - 1 ]; 
} /* end ‘if found’ */ 
else 
pDistanceQ = ( *pDistanceQ ).pNext; 
}; 
}; 


(*CV[ c 1 ).days = (*CVI_c] ).days + (*CV[ c] ).stops; 
(*CV[c] ).cP = 1; 
(*CV[E.c] ).RV = 1; 

); /* end CV loop*/ 


} /*** end of procedure DetermineDistance ***/ 


/“e Procedures and Functions for Heap Sort *********/ 


short Left( short i ) 


/* 
Called from procedure heapify. Determines index 
for left node. 
*/ 
{ 
return( 2 * i ); 
} /*** end of procedure Left ***/ 


short Right( short i ) 


/* 
Called from procedure heapify. Determines index 
for right node. 
Bs i 
{ 
return( (2*i)+1); 
) /*** end of procedure Right ***/ 


short Exchange( short i, short largest ) 


/* 
Called from procedure Heapify. Exchanges the 
nodes in the array. 
*/ 


struct ship “temp; 
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temp = CV[i J; 
CV[ i] = CV[ largest J; 
CV[ largest ] = temp; 


) /*** end of procedure Exchange ***/ 


void Heapify( short i, short heapSize ) 
/* 


Called from procedure SortQueue, and recursively 
te astablish heap array. 
*/ 


short 1, r, largest; 


1 = Left i ); 

r = Right( i ); 

if( (1 <= heapSize ) && ((*CV[1] ).status > (*CV[i] ).status ) ) 
largest = 1; 

else 
largest = i; 


if( (r <= heapSize ) && ( ( *CV[ r] ).status > ( *CV[ largest ] ).status ) ) 
largest = r; 


if( largest != i ) 
{ 
Exchange( i, largest ); 
Heapity( largest, heapSize ); 
}; 


} /*** end of procedure Heapify ***/ 


void SortQueue() 
/* 
Called from main. Organizes Dial solution into heap array, then 
sorts the array by status first, using a heap sort, then by 
distance by using linear comparisons. 


*/ 
short, i, heapSize; 


heapSize = numCV; 
c = heapSize/2; 


while( c >= 1) 

Heapify( c - 1, heapSize - 1 ); 
}; Ee 

c = heapSize; 
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while( c >= 2 ) 
{ 


Exchange( 0, ¢ - 1); 
heapSize--; 

Heapify( 0, heapSize - 1 ); 
c-; 


i 


for( c= 0; ¢ <= numCV - 2; c++ ) 
if( ( *CV[ c] ).status == (*CV[ c+ 1] ).status ) 
if( (*CV[_c ] ).distance > (*CV[_c + 1 ] ).distance ) 


{ 
Exchange(c,c + 1); 
} 
else 
if( ( *CVI c} ).distance == (*CV[ c + 1 ] ).distance ) 
Exchange(c,c +1); 
}; 
shipQ = CVI 01; 
for( c = 0; c <= numCV - 2; c++ ) 
{ 
(*CVI_c] ).pNext = CVI c+ 1); 
(*CV[.c +1] ).pBack = CVI c); 
(*CV[ 0] ).pBack = NULL; 
) /*** end procedure SortQueue ***/ 


jews" Procedures and Functions for display ******/ 


void ChangeQueue() 
Called to change the order of the ship queue. 
*/ 


short ID, pos, c = 0; 
struct ship “*pTemp; 
Boolean good = false; 


while( !good ) 
{ 
printf( "\nWhich ship would you like to change? “ ); 
scanf( "%d", &ID ); 
if((ID<0) 11 (ID>numCV )) 
printf( "Not a valid entry. Please re-enter. \n" ); 
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void 


\ 


good = false; 
while( !good ) 
{ 


printf( "\nWhere should I put it in the queue? " ); 
scanf( "%d", &pos ); 
if({ pos <0) 11 (pos > numCV )) 
printf( "Not a valid entry. Please re-enter. \n” ); 


good = true; 


else 
}; 


while( (*CV[ c} ).ID != ID ) 
C++; 

(*CV[ c- 1] ).pNext = (*CVI c] ).pNext; 
(*CV[_c + 1] ).pBack = (*CV[c ] ).pBack; 
(*CVI pos - 2] ).pNext = CVI c ]; 

(*CV[ pos - 1 J ).pBack = CV[c ]; 

(*CV[ c] ).pNext = CV[ pos - 1]; 

(*CVI c ] ).pBack = CVI pos - 2); 


if( pos = 1) 
shipQ =CVic J; 
pTemp = shipQ; 
for( c = 0; c <= numCV - 1; c++) 
{ 
CV[ cj = pTemp; 
pTemp = ( *pTemp ).pNext; 
}; 


/*** end of procedure ChangeQueue ***/ 


ChangeStatus() 


short ID, state, c = 0; 
Boolean good = false; 


while( !good ) 
{ 


printf( "\nWhich ship would you like to change? " ); 
scanf( "%d", &ID ); 
if((ID <0) 1! (ID>numCV )) 
printf( "Not a valid entry. Please re-enter. \n" ); 
else 
good = true; 
; 


good = false; 
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void 


void 


while( !good ) 
{ 


printf( "\nWhat is its new status? " ); 
scanf( "%d", &state ); 
if( (state< 1) 11 (state>5)) 
printf( "Not a valid entry. Please re-enter. \n" ); 


good = true; 


else 
}; 


while( (*CV[ c } ).ID != ID * 
C++; 
if( (*CV[ c ] ).status == 4 ) 
(*CV[ c] ).days = ((*CV[c]).distance / 24 / xsitSpeed) + ("CV[c]).stops; 
(*CV[_c ] ).status = state; 
SortQueue(); 


/*** end of procedure ChangeStatus ***/ 


ChangeStops() 
short ID, stops, c = 0; 
Boolean good = false; 


while( !good ) 
{ 
printf( "\nWhich ship would you like to change? ” ); 
scanf( "%d", &ID ); 
if(( ID <0) Il (ID>numCv )) 
printf( "Not a valid entry. Please re-enter. \n" ); 


good = true; 


else 
); 


printf( "\nHow many days in port? " ); 
scanf( "%d", &stops ); 
while( (*CV{c]).ID != ID ) 
C++; 
(*CV[ c] ).days = (*CV{ cj ).days - (*CV[c ] ).stops; 
("CVI c] ).stops = stops; 
(*CV[c] ).days = (*CV[ cc] ).days + (*CV[c ] ).stops; 


/*** end of procedure ChangeStops ***/ 


ChangeMRT() 
/* 
Called from procdure GenerateTable. Allow user 
to change unit make-ready times. 
*/ 
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void 


void 


short ID, mrt, c = 0, steamingDays; 
Boolean good = false; 
while( {good ) 
printf( "\nWhich ship would you like to change? " ); 
scanf( "%d", &ID ); 


if((ID <0) 11 (ID >numCV )) 
printf( "Not a valid entry. Please re-enter. \n" ); 


good = true; 


else 
}; 


printf( "\nWhat is the new make-ready time (in days)? " ); 
scanf( "%d", &mrt ); 


while( ( *CV[ c] ).ID != ID ) 
C++; 

(*CV[ c J ).days = mrt; 

steamingDays = ((( *CV[ c ] ).distance )/xsitSpeed } / 24; 

if( ( (*CV[_c] ).status == 4 ) && ( ( *CVI c] ).location != source ) ) 
(*CV[ c} ).days = 30; 

else 
(*CV[ c] ).days = (*CV[ c ] ).days + steamingDays; 


/*** end of procedure ChangeMRT ***/ 


ChangeSpd() 
/* 
Called from procedure GenerateTable. Allows user 
to change the standard response speed of 18 kts. 
*/ 


short c, spd, steamingDays; 


printf( “\nWhat is the new transit speed (in knots)? " ); 

scanf( "%d", &spd ); 

for( c= 0; c <= numCV - 1; c++) 

( 
steamingDays = ( ( ("CVI c ] ).distance )/xsitSpeed ) / 24; 
(*CV[_c] ).days = (*CV[ c ] ).days - steamingDays; 
steaming Days = (( (*CV[c } ).distance )/spd ) / 24; 
(*CV[ c ] ).days = ( *CV[ c ] ).days + steamingDays; 

}; 


/*** end of procedure ChangeSpd ***/ 


GenerateTable() 


48 


Called from main. Generates a table representing the 
ship queue in tabular form, in order of preference. 
Allows user to modify parameters as necessary. 


st 


short c, opt; 
char chg; 
Boolean valid = false; 


good = false; 
while( !good ) /* while not happy with the table... */ 
{ 


printf( "CV \t Status \t Location \t Distance \t Port Days \t Days \n" ); 
for( c = 0;c <= numCV - 1; c++ ) 
{ 
printf("%d \t %d \t\t %d \t\t"(°CVIc]).ID,(CCVIc]}).status, 
(*CV[c]).location); 
printf('" Yd \t\t %d \t\t %d \n", (*CV[_c ] ).distance, 
(*CV[ c] ).stops, (*CV[ ¢ ] ).days ); 


valid = false; 
while( !valid ) /* keyboard input isn't valid! */ 
{ 
printf( '"\nWould you like to make a change? " ); 
FlushFile( stdin ); 
scanf( "%c", &chg ); 
valid = ( (chg == ‘Y‘) {1 (chg ==’y’) 1 (chg == 'n’) || (chg == 'N’) ); 
if( !valid ) 
printf( "\nNot a valid entry. Please try again. \n" ); 
if( (chg =="Y') I 1 (chg =='y')) 
{ 


printf( ""\n\t{ 1 }\tChange Ship Queue\n" ); 
printf( "\tl 2 ]\tChange Ship Status\n" ); 
printf( ""\t{ 3 ]\tChange Number of Port Stops\n" ); 
printf( “\t{ 4 ]\tChange Ship Make-ready Time\n” ); 
printf("\tl 5 ]\tChange Transit Speed\n" ); 
printf( "\tl 6 ]\tNo Change\n" ); 
printf( "\nSelect option: “ ); 
scanf( "%d", &opt ); 
switch( opt ) 
{ 
case 1: 
ChangeQueue(); 
break; 


case 2: 
ChangeStatus(); 
break; 


case 3: 
ChangeStops(); 


49 


break; 


case 4; 
ChangeMRT(); 
break; 
case 5: 
ChangeSpd(); 
break; 
default: 
good = true; 
); /* end of switch case */ 
} /* end of if change = yes */ 
else 
good = true; 
); /* end of change loop */ 
}; /* end ‘while good' loop */ 
} /** end of procedure GenerateTable ***/ 


/**** Procedures for calculating MOEs ******/ 


short FindDelay( short cv ) 


/* 
Called from procedure CalculateMOE. Determines 
delay, in days, for FDO to arrive in the region. 
This value is returned to CalculateMOE. 
*/ 
{ 
short t, delay = 0; . 


for( t = 0; t <= cv - 1; t++ ) 
{ 
if( delay <= ( *CV[ t ] ).days ) 
delay = ( *CV[ t ] ).days; 
}; 


return( delay ); 
) /*** end of procedure FindDelay ***/ 
void OutputMOE2( short ships, short delay, FILE “fp ) 
/* 


Called from procedure CalculateMOEs. Calculates the 
total percentage of the FDO requirements on station 
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void 


each day. This data is output to an external ASCII 
text file called ECRM Output and can be opened by 
commercially available software for graphical analysis 


of the response. 
*/ 
short d=0,t; 
float total = 0; 
while( d <= delay ) 
{ 
for( t = 0; t <= ships - 1; t++ ) 
{ 
if( (*CV[ t] ).days == d ) 
total = total + (*CVI t ] ).cP; 
}; 
fprintf( fp, "%f\t", total/ships ); /* total cP or hV required */ 
d++; 
}; f? while the day counter is less than the total delay for the 


FDO to arrive. */ 
fprintf( fp, "\n\n\n" ); 


/** end of procedure OutputMOE2 ***/ 


CalculateMOEs() 
/* 

Called from main. Calculates the total number of days 
required for the total FDO package to arrive. This 
information is displayed on screen. The total percentage 
of the FDO package on station each day is output to 
an external file. Also provides data on the FDO package 
one level below, and one level above the recommended 
levels. 


*/ 


short delay, cv, cvLow, cvHigh, d, m; 
FILE ‘*fileOut; 


fileOut = fopen( "ECRM Output”, "w" ); 


cv = div( fDO, 10 ).quot; 
if(cv == 1) 

cvLow = cv; 
else 

cvLow = cv - 1; 
cvHigh = cv + 1; 
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printf( "\nThis crisis requires %d carrier", cv ); 
if(cv == 1) 
printf( ".\n" ); 
else 
printf( “s.\n" ); 
printf( "Delay until complete FDO arrives is %d days. \n", delay = FindDelay( cv ) ); 
OutputMOE2( cv, delay, fileOut ); 


printf( "For %d carrier", cvLow ); 
if( cvLow == 1) 
printf(", " ); 
else 
printf( "s, " ); 
printf( "the delay is %d days.\n", delay = FindDelay( cvLow ) ); 
OutputMOE2( cvLow, delay, fileOut ); 


printf( "For %d carriers, the delay is %d days.\n", cvHigh, delay = FindDelay( 
cvHigh ) ); 

OutputMOE2( cvHigh, delay, fileOut ); 

fclose( fileOut ); 


} /*** end of procedure CalculateMOEs ***/ 
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[TEASER EE EERE ES ROR EEE OSES SENSE ESET NEE OH EEE EENE SEAN EE EES EHS OH OSE 8: 


Dial's Algorithm Module : 
Written by LT Takashi R. Yamamoto * 
Operations Research Department, Naval Postgraduate School * 


* 


Called from main to implement the shortest-path . 


* 
* 
* 
* 
* 
* * 


TERR ER AREER ERA RER TESS SORENESS TEER ESR AERTS ARON EEE EERE EEE SAS OEE / 


#include <limits.h> 
#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include "Globals.c" 


J setasananeenesnee PROCEDURES AND FUNCTIONS SRUSTENSR RATES EEE / 


void BuildAdjList( struct edge *pNew ) 
/* 
Called from procedure ReadNetworkData to build 
and update a sparse matrix composed of nodes 
representing adjacent geographic regions. The 
adjacent nodes are assumed to be normal 


transit legs taken by ships. 
{ 
if( adjList[ ( *pNew ).nodeOut - 1 ].pAdjNodes == NULL ) 
adjList[ ( *pNew ).nodeOut - 1 ].pAdjNodes = pNew; 
else 
{ 
( *pNew ).pNext = adjList[ ( *pNew ).nodeOut - 1 ].pAdjNodes; 
adjList{ (*pNew ).nodeOut - 1 ].pAdjNodes = pNew; 
}; 
) /™* end BuildAdjList ***/ 


void MakeNodeSet( struct nodeSet *pVert ) 
/* 
Called from procedure ReadNetworkData to build 
a linked list of nodes representing geographic 
regions in the Pacific theater. 


a 
{ 
if( pFirst == NULL) /* if it’s the first node */ 
{ 
pFirst = pVert; 
pLast = pFirst; 
} 
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void 


void 


else /* all other nodes that follow */ 
( *pLast ).pNext = pVert; 
(*pVert).pBack = pLast; 
pLast = ( *pLast ).pNext; 

); 


/*** end MakeNodeSet ***/ 


CreateBucket( short source ) 


/* 
Called from procedure ReadNetworkData to build 
an array of "buckets" where the index represents 
the distance from the source node to other regions. 
of 
struct nodeSet “pIndex, *pIndex2; 


pindex = pFirst; 
if( ( *pIndex ).node == source) /* if first node is source node */ 


B[ maxEdge ].pTr = ( *pIndex ).pNext; 
( *( *pindex ).pNext).pBack = NULL; 


else /* source node is not first node */ 


B[ maxEdge ].pTr = pIndex; 
while( ( *pIndex ).node != source ) 
pindex = ( *pIndex ).pNext; 
(*( *pIndex ).pBack ).pNext = ( *pIndex ).pNext; 
( *( *pIndex ).pNext ).pBack = ( *pIndex ).pBack; 
pindex2 = ( *pIndex ).pBack; 
); 


( *pIndex ).pNext = NULL; 

( *pIndex ).pBack = NULL; 

B[ 0 ].pTr = pIndex; /* places source region into */ 

(*B[ 0 ].pTr ).distance = 0; /* first location. */ 


( *B[ 0 J.pTr ).bucket = 0; 
if( ( *pIndex2 ).pNext == NULL ) 
pLast = pIndex?2; 


/*** end CreateBucket ***/ 


ReadNetworkData( short source ) 


/* 


struct 


short 
struct 


while( ! 


}; 


Called from procedure Dial to read in the 
graph data representing regions in the 
Pacific, and the DMA distances associated 
with traveling between those points. 


“4 


nodeOut, nodeln, i, dist; 
edge “*pTemp; 


(feof( data ) ) ) 


fscanf( data, "%d %d %d", &nodeOut, &nodeln, &dist ); 
pTemp = malloc( sizeof( struct edge ) ); 
(*pTemp).nodeOut = nodeOut; 
(*pTemp).nodelIn = nodeln; 
(*pTemp).distance = dist; 
(*pTemp).pNext = NULL; 
BuildAdjList( pTemp ); 

pTemp = malloc( sizeof( struct edge ) ); 
(*pTemp).nodeOut = nodelIn; 
(*pTemp).nodelIn = nodeOut; 
(*pTemp).distance = dist; 
(*pTemp).pNext = NULL; 

Build AdjList( pTemp ); 


for(i = 1; i <= maxNodes; i++ ) 


{ 


); 


pVertex = malloc( sizeof( struct nodeSet ) ); 
( *pVertex ).node = i; 

( *pVertex ).distance = infinity; 

( *pVertex ).pNext = NULL; 

( *pVertex ).pBack = NULL; 

( *pVertex ).bucket = maxEdge; 

P{ i] = pVertex; 

MakeNodeSet( pVertex ); 


CreateBucket( source ); 
fclose( data ); 


/*** end of ReadData ***/ 


/*** procs and funcs for Dial’s ***/ 


nodeSet 


*GetNode( short i ) 


/* 
Called from procedure FindMiIn to “pull” the 
next node from the distance bucket. This node 
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represents the next closest region. 
*/ 

struct nodeSet *pTemp; 

pTemp = B[ i ).pTr; 

Bl i |.pTr = ( *B[ i }.pTr ).pNext; 

return( pTemp ); 


} /** end of GetNode ***/ 


void Update( struct nodeSet *V ) 
/* 


Called from procedure FindMin to update the 
distance bucket and pointer array after getting 


the next node. 
*/ 
( 
struct edge ‘J; 
struct nodeSet *marker; 
short mod; 


J = adjList[ (*V).node - 1 ].pAdjNodes; 
while( J != NULL ) 
{ 


if( ( (*P[ (*J).nodeln ]).distance ) > 
( (*V).distance + (*J).distance ) ) 
{ 


/*If distance in P is greater than V.distance + J.distance */ 
(*P[ (*J).nodeln ]).distance = ( (*V).distance + (*J).distance ); 
mod = fmod( (*P[ (*J).nodelIn ]).distance, maxEdge ); 

marker = B[ mod ].pTr; 


if( B[ mod ].pTr == NULL ) /* Bucket is empty */ 
{ 


if( *P[ (*J).nodeln ]).pBack == NULL ) 


BE (*P[ (*J).nodeln }).bucket ).pTr 

= (*P[ (*J).nodeln }).pNext; 
(*(*P[ (*J).nodeln }).pNext ).pBack = NULL; 
(*P[ (*J).nodelIn ]).pNext = NULL; 
BI mod ).pTr = &(*P[ (*J).nodeln J); 
(*P[ (*J).nodeln }).bucket = mod; 


else /* not first node in line */ 


( *(*PL (*J).nodeln }).pBack ).pNext 


= (*P[ (*J).nodeln )).pNext; 
( *(*P[ (*J).nodeln }).pNext ).pBack 

= (*P[ (*J).nodeln ]).pBack; 
(*P[ (*J).nodeln }).pNext = NULL; 
(*P[ (*J).nodeIn }).pBack = NULL; 
B[ mod J.pTr = &(*P[ (*J).nodein }); 
(*P[ (*J).nodeln ]).bucket = mod; 


) /* end of "if bucket empty” */ 
else /* Bucket not empty */ 
{ 


while( ( *marker ).pNext f= NULL ) 
marker = ({ *marker ).pNext; 


if( (*P[ (*J).nodeln }).pBack == NULL ) 

{ 
B[ (PI (*J).nodeln }).bucket ].pTr 

= (*P[ (*J).nodeln }).pNext; 

( *CP[ (*J).nodeIn ]).pNext ).pBack = NULL; 
(*P[ (*J).nodein ]).pNext = NULL; 
( *marker ).pNext = &(*P[ (*J).nodeln )); 
(*P{ (*J).nodeln ]).pBack = marker; 
(*P[ (*J).nodeln }).bucket = mod; 


} 
else /* not first node in line */ 
{ 
( *(*P[ (*J).nodein J).pBack ).pNext 
= (*P[ (*J).nodeln }).pNext; 
( *(*P[ (*J).nodeln }).pNext ).pBack 
= (*P[ (*J).nodeln }).pBack; 
(PL (*J).nodeln }).pNext = NULL; 
(*P{ (*J).nodein J).pBack = NULL; 
( *marker ).pNext = &(*P[ (*J).nodeln J); 
(*P[ (*J).nodeln ]).pBack = marker; 
(*P{ (*}).nodein }).bucket = mod; 
}; 
}; /* end of "if bucket is empty */ 
}; /* end of "if P[d] > V[d] + J{d]" */ 
J = ("J ).pNext; 
; /* end of "while" loop */ 


) /*** end Update ***/ 


void FindMin(Q 
* 


Called from procedure Dial. Calls function GetNode 
to find next nearest node. Then updates the set S 


void 


which is a linked list of the nodes in order of 
distance from the source. This list is also stored 
in the pointer array P. Lastly, calls the procedure 
Update to update the distance bucket and pointer 
array after obtaining a node. 


af 


short x; 
struct nodeSet "Vertex; 


while( B[ (x = fmod( D, maxEdge ) ) |.pTr == NULL ) 
D++; /* looking for a bucket that isn't empty */ 

Vertex = GetNodef x ); 

if( setS == NULL ) 

{ /* updates set S */ 
setS = Vertex; 
sLast = setS; 


else 


( *sLast ).pNext = Vertex; 
sLast = ( *sLast ).pNext; 
); 
Update( Vertex ); 


/*** end of FindMin ***/ 


Dial() 
/* 
Called from ECRM main to implement Dial’s 
Algorithm, a modification of the Dijkstra 
algorithm, to obtain the shortest-path 
distances from a source node to all nodes 
in the network. 
at 
short count; 
data = fopen( "Network", "r" ); 
ReadNetworkData( source ); 
D=0; 


count = 0; 
while( count != maxNodes ) 


count++; 
FindMin(); 
); 


/*** end Dial ***/ 
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¢ USN Surface Ship Employment and Crisis Response Model . 
* * 
Gs Written by LT Takashi R. Yamamoto . 
* Operations Research Department, Naval Postgraduate School * 
. March, 1993 x 
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#include <limits.h> 
#include <stdio.h> 
#include <stdlib.h> 
#include "structs.h" 


[Press treenseeres GLOBAL VARIABLES AOESSTAPERR TEESE S: / 


FILE “data; 

short crisis, region, source, D, i, xsitSpeed; 
short fDO; 

Boolean good; 

struct node  adjList{ maxNodes ]; 

struct pNode B[ maxEdge + 1]; 


struct nodeSet *P[ maxNodes + 1 ]; 
struct nodeSet *pFirst, *pLast, *setS, *sLast, *sMarker; 
struct nodeSet *pVertex; 


struct Posit CVLocation{ numCV ]; 
struct ship “*CV[{numCV]J; 


f[ PARSER ANNE NST R aS MAIN CODE TEST AN atone rene eeeeeeee / 


main() 
{ 

InitializeQ; 

ReadLocation(); /* Reads in CV location matrix and stores 
the data in a linked list. +] 

DetermineCrisis(); f Gets crisis code and region code 
Currently limited to a single crisis 
and region. */ 

DialQ; <* Calls Dial's Algorithm to compute the 
shortest-path distances from the source 
region tc other regions in the theater. */ 

DetermineDistance(); ;* Establishes distance queue, updates ship 
type data structures for each ship. */ 

SortQueue(); * Creates ship queue, then uses heap sort to 


sort queue by status and distance. */ 


GenerateTable(); /* Generates output of queue options on 
console screen. User can modify response 
* 


queue as necessary. / 


CalculateMOEs(); /* Calculates measures of effectiveness. */ 


an OF Ff WY NY 


on & WO NO —& 


APPENDIX B 


GEOGRAPHIC LOCATION CODES 


Eastern Pacific ( San Francisco, San Diego, Seattle ) 
Mid-Pacific ( Pearl Harbor ) 

Northwestern Pacific ( Yokosuka ) 

Southwestern Pacific ( Philippines ) 

Southwestern Pacific ( Singapore ) 


Indian Ocean 


SHIP STATUS CODES 


In Transit / Local Operations 
Upkeep 

Returning from Deployment 
Not Available ( Employed ) 


Major Maintenance Availability 
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APPENDIX C 


HYPOTHETICAL FDO CRISIS CODES 


Crisis Cod Tit Descripti FDO Cod 


1 Coup d’Etat Coup in “friendly” nation 10 
2 Bombing Bombing of embassy /base 10 
3 Revolution Revolution in country 10 


of interest 


4 Civil War Civil War in country 10 
of interest 


5 Territorial Dispute Friendly nation involved 20 
LRC/LIC 
6 Hostages US/ Allied hostages 21 
7 Attack on US Forces MRC 32 
8 Humanitarian Aid Non-disaster various 
9 Non-combatant Rescue of citizens 11 
Evacuation 
10 Blockade UN Maritime Interdiction 30 
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